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We study by means of the recursive Green's function technique the local density-of-states of (finite 
and semi-infinite) multi-band spin-orbit coupled semiconducting nanowires in proximity to an s-wave 
superconductor and attached to normal- metal electrodes. When the nanowire is coupled to a normal 
electrode, the zero-energy peak, corresponding to the Majorana state in the topological phase, 
broadens with increasing transmission between the wire and the leads, eventually disappearing for 
ideal interfaces. Interestingly, for a finite transmission a peak is present also in the normal electrode, 
even though it has a smaller amplitude and broadens more rapidly with the strength of the coupling. 
Unpaired Majorana states can survive close to a topological phase transition even when the number 
of open channels (defined in the absence of superconductivity) is even. We finally study the Andreev- 
bound- state spectrum in superconductor- normal metal-superconductor junctions and find that in 
multi-band nanowires the distinction between topologically trivial and non-trivial systems based on 
the number of zero-energy crossings is preserved. 

PACS numbers: 74.78.Na,74.45.+c,71.10.Pm,74.78.-w 



I. INTRODUCTION 

Since the first prediction^ of real solutions to the Dirac 
equation, known as Majorana fermions, there have been 
many attempts to demonstrate their occurrence in na- 
ture^, but a clear evidence is still lacking. Besides 
the natural search for these elusive particles in high- 
energy physics, it has been recently suggested that Ma- 
jorana fermions can exist as exotic excitations in certain 
condensed-matter systems^'"^. Such solid-state realiza- 
tions include fractional quantum Hall states at filling fac- 
tor V = 5/2^, p-wave superconductors and superfluids^'^, 
three-dimensional topological insulators in proximity to 
5- wave superconductors^, as well as spin-orbit cou- 
pled semiconductors in a magnetic field with proximity- 
induced 5-wave superconducting pairing^^^. The impor- 
tance of finding Majorana fermions in condensed-matter 
systems is not only related to their fundamental interests. 
It is also rooted in the non-Abelian braiding statistics of 
these particles which could be exploited as a basis for 
decoherence-free topological quantum computation^^. 

In this paper we focus on a specific proposal real- 
ized with spin-orbit-coupled semiconducting nanowires 
in proximity to an s-wave superconductor (S) and sub- 
jected to an in-plane magnetic field ^^^^ (a system that, 
for the sake of simplicity, will be henceforth termed "S- 
nanowire"). This wire can support Majorana-fermion 
bound states at its ends when parameters such as chem- 
ical potential, magnetic field, and superconducting pair- 
ing are properly tuned^^. The S-nanowire is said to be in 
the topological phase when Majorana bound states are 
present, while it is topologically trivial otherwise. 

In order to assess the presence of Majorana fermions in 
such solid-state systems, it is of primary importance to 
predict clear signatures of the topological phase, which 
could be then used to guide experiments. A very relevant 
quantity is the Local Density-of-States (LDOS), which 



can be accessed in scanning tunneling microscopy. LDOS 
calculations have already been carried out in Refs. 17 
and 18, but were restricted to finite-size S-nanowires in 
which superconductivity is uniform along the wire. The 
influence of normal segments inside a one-dimensional 
S-nanowire was investigated in Ref. 19. In this work 
we exploit the recursive Green's function technique to 
study the LDOS of both finite and semi-infinite multi- 
band S-nanowires in the presence of normal electrodes 
(NS junctions), as well as Supercoductor-Normal metal- 
Superconductor (SNS) junctions. In this manner we have 
been able to i) reproduce previous results^ ii) gen- 
eralize them to account for the coupling to normal elec- 
trodes and/or for a finite width of the wire, and iii) con- 
sider more complicated structures like SNS junctions. 



The paper is organized as follows. In Sec. II we de- 
scribe the model and we introduce the corresponding 
Hamiltonian. The phase diagram is presented and new 
analytical expressions for the phase boundaries are re- 
ported. Section III is devoted to the discussion of our 
numerical results for the LDOS. We start in Sec. Ill A 
with the case of a finite-size multi-band S-nanowire^^. 
We then consider a single NS junction between semi- 
infinite leads in Sec. IIIB. This situation can arise for 
instance when the nanowire is only partially in proxim- 
ity to a bulk superconductor so that part of the nanowire 
is in the normal state. Here we also study the evolution 
of the LDOS across the topologically trivial/topologically 
non-trivial phase transition. In Sec. Ill C we turn our at- 
tention to a SNS junction, illustrating numerical results 
for the LDOS and for the Andreev-bound-state spectrum. 
Finally, in Section IV we draw our main conclusions. 
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II. MODEL HAMILTONIAN 

Our calculation of the LDOS is based on a recursive 
Green's function method especially designed for tight- 
binding Hamiltonians. In order to apply this method 
to the present case, we describe the semiconducting 
nanowire as a square lattice with a finite width W in 
the direction perpendicular to the axis of the wire (the 
^-direction) and lattice constant a. Assuming the pres- 
ence of Rashba-type spin-orbit (SO) couphng, of strength 
a, and a Zeeman field V along the wire (;r-direction), 
the lattice discretization of the usual continuum-model 
Hamiltonian^^'^^'^^ reads^^ 

(i,j),cr i,cr 

+v J2 <^%A,A,a' ■ (1) 

Here £q = 4t is a uniform on-site energy which sets the 
zero of energy, a* are spin- 1/2 Pauli matrices, Vij = x • 
dij, and u^j = y • dij, with dij = (r^ -rj)/\ri- rj \ being 
the unit vector connecting site j to site i. 

If we now allow sections of the nanowire to be in con- 
tact with a bulk 5-wave superconductor, the proximity 
effect induces a non-vanishing superconducting pairing 
in these sections so that the complete Hamiltonian be- 
comes 



where 



'H = 'Hrz + 'Hs 



(2) 



(3) 



For simplicity we assume A{i) to be piecewise constant, 
with |A(z)| = A in the regions in contact with the super- 
conductor and A{i) = otherwise. Moreover, we assume 
that a barrier is present at the boundary between proxi- 
mized and non-proximized sections, leading to a decrease 
of the value of the hopping energy t and of the SO cou- 
pling a by the same factor 7.^^ 

Finally, it is convenient to introduce the Nambu 
spinors = (q^^, q^^, ^, — ^)-^ and re-write the 
Hamiltonian (2) in the form 



where 



A{i)S,j -ay h^^^iij) a'^ 



(4) 



(5) 
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FIG. 1: (Color online) Phase diagram of an infinite supercon- 
ducting wire as a function of the chemical potential /x and the 
Zeeman field V. White regions correspond to a trivial sys- 
tem (Q = +1), while dark regions identify the topologically 
non-trivial phase (Q = ~1)- R^d thick lines illustrate the 
prediction in Eq. (8) for the phase boundaries. These results 
refer to the following set of parameters: W/a = 10, a/t = 0.1, 
and A/t = 0.1. 



is the Bogoliubov-de Gennes (BdG) Hamiltonian^^. 

The tight-binding Hamiltonian (5) will be the starting 
point of our analysis. 



A. Phase diagram 

Before proceeding with the study of the LDOS we need 
to know in which regions of parameter space we should 
expect Major ana fermions. This problem has been ad- 
dressed in Refs. 13 and 14 in the strictly one-dimensional 
case and in the continuum limit (a ^ 0): the S-nanowire 
is in the topological phase when > while it 

is in the trivial phase otherwise. Thus the phase bound- 
ary occurs along the line implicitly defined by 



(6) 



The phase diagram of multiband {W/a ^ 1) nanowires 
has been investigated numerically in Refs. 17 and 20. The 
phase of a uniform system can be determined for instance 
from the evaluation of the following Pfaffian formula^^ for 
the topological invariant 

Q = sign{Pf[HBdG(0)aV^] Pf [HBdG(V«)aV^]} , (7) 

where HBdcikx) is the Fourier transform of the BdG 
Hamiltonian in Eq. (5), while are Pauli matrices act- 
ing on the particle-hole degrees-of-freedom. In Fig. 1 we 
report numerical results for a given system {W/a = 10, 
a/t = 0.1, and A/t = 0.1) as a function of the chemical 
potential /i and Zeeman field obtained using an algo- 
rithm developed by Wimmer^^. The topologically trivial 
phase corresponds to Q = +1 (white regions), while the 
non-trivial one is signaled by Q = — 1 (dark regions). 
The phase boundaries in this figure (red thick lines) have 
been derived analytically and are given by the following 
result: 



{fi-eo-ex± 2tf + A' 



(8) 
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FIG. 2: Experimental setup adopted in measurements of the 
LDOS of a nanowire (in this case comprising a normal elec- 
trode and a superconducting segment separated by a bar- 
rier). The differential conductance obtained from the current 
between the tip of the scanning tunneling microscope (STM) 
and the nanowire is proportional to the LDOS of the nanowire 
according to Eq. (10). 



Here £0 + T 2t are the eigenenergies of HBdcik) for 
k = 0,7r/a, respectively, when V = A = fi = 0. The 
following expression holds for the energies ^a- 



Att 



A = 1, . . . , n = W/a . 

(9) 

A thorough derivation of Eqs. (8) and (9) is given in 
Appendix A. 
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FIG. 3: (Color online) (a) LDOS of an isolated superconduct- 
ing nanowire in the topologically non-trivial phase (/i/t = 0, 
V/t = 0.2) at an energy very lose to the chemical potential 
(^ ^ 0). Bound states at both ends of the wire are apparent, 
(b) LDOS at a given position (x/a = 4, ^//a = 5) as a function 
of energy. A sharp peak corresponding to a Majorana bound 
state is present at ^ = 0. (c) LDOS at ^ as a function 
of X along the middle of the wire {y/a = 5). 



III. LDOS - NUMERICAL RESULTS 

As we mentioned in the Introduction, the LDOS 
can be accessed in experiments using scanning tunnel- 
ing microscopy (STM). The STM experimental setup is 
sketched in Fig. 2. When the metallic tip of the STM 
is moved close to the nanowire a tunneling current can 
flow. By locally measuring this current / as a function 
of the tip-sample bias voltage V, the LDOS at a given 
position r and energy E can be reconstructed from the 
differential conductance at eV = 

^{r,eV)(xJV{r,E = eV) . (10) 

As a consequence, it is particularly interesting to inves- 
tigate the LDOS theoretically and make predictions that 
can be tested in experiments. 

In what follows the LDOS is computed through the 
standard relation 

Af{r,E) = -^3m {Tr[G{r,E)]} , (11) 

where G(r, E) is the Green's function and the factor 2 in 
the denominator is introduced to avoid a double count- 
ing of particle and hole degrees-of-freedom intrinsic in 
the BdG formalism. We have computed G{r,E) using 
a recursive Green's function technique similar to the one 



adopted in Ref. 17, suitably generalized to include the ef- 
fects of semi-infinite leads^^. For simplicity, in the follow- 
ing we fix the width W = 10 a, the SO coupling strength 
a = 0.1 t and the superconducting pairing A = 0.1 t. 



A. Isolated S-nanowire 

Let us first consider an isolated S-nanowire of finite 
length (L = 100 a). This situation has been already 
addressed before^^'^^ and it is considered here for the 
sake of reference. We consider two cases: i) /i = and 
V/t = 0.2 (with one open channel in the absence of su- 
perconducting pairing) and ii) /i = and V/t = 0.6 (with 
two open channels). According to Fig. 1, in case i) the 
wire is topologically non-trivial, while in case ii) the wire 
is topologically trivial. For case i). Fig. 3(a) shows that 
the LDOS at an energy close to the chemical potential 
is characterized by the presence of bound states at both 
ends of the wire. The presence of Majorana bound states 
appears as a sharp peak at zero energy in the LDOS 
as a function of energy at a given position in space [see 
Fig. 3(b)]. According to Fig. 3(c), these Majorana bound 
states have oscillating wavefunctions which decay expo- 
nentially inside the bulk of the S-nanowire with a typical 
length scale (effective superconducting coherence length) 
^ ^ 10 a. 

On the contrary, in case ii) where the wire is topo- 
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FIG. 4: (Color online) (a) LDOS of an isolated supercon- 
ducting nanowire in the topologically trivial phase (/x/t = 0, 
V/t = 0.6) at E/t 0.002. (b) LDOS at a given position 
(x/a — 2, y/a — 5) as a function of energy. Owing to the 
coupling between the two Majorana end states, two Dirac- 
fermion modes appear at finite energy [E/t ^ 0.002). (c) 
LDOS at E/t ^ 0.002 as a function of x along the middle of 
the wire (y/a = 5). 
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FIG. 5: LDOS of a superconducting nanowire coupled to 
a normal lead as a function of energy, at fixed positions in 
space, in the non-trivial [/n/t = and V/t = 0.2, panel (a) 
and (c)], and trivial [/i/t = 0, V/t = 0.6, panel (b) and (d)] 
phase. Panels (a) and (b) refer to a position just inside the 
superconducting part (x/a — 51), while panels (c) and (d) 
refer to a position just inside the normal part of the junction 
(x/a — 50). The interface is at x/a = 50.5. 



logically trivial (and presents two transverse channels in 
the absence of superconducting pairing) the LDOS at the 
chemical-potential energy is almost zero throughout the 
wire, while it shows spatial features only at finite ener- 
gies [see Fig. 4(a) for E/t ±0.002, where E is measured 
from the chemical potential]. In particular, as shown in 
Fig. 4(c), the LDOS oscillates and decreases moving to- 
ward the center of the wire, the length scale of the expo- 
nential drop, ^ ~ 30 a, being much larger than for case i). 
Now, if the two channels were decoupled, two Majorana 
modes would have appeared at each end of the wire, one 
for each open channel. However, since the two transverse 
channels are actually coupled in the wire, a single fermion 
at each end appears at a finite energy as though com- 
ing from the hybridization of the two Majorana modes. 
Fig. 4(b) shows one pair of peaks at E/t ±0.002 and 
another pair (with smaller amplitude) at E/t ±0.007. 
The presence of two pairs of peaks is due to the long co- 
herence length which allows the two localized fermions 
at the ends of the wire to strongly hybridize, lifting the 
parity degeneracy of the system. We have checked that 
for a longer wire (L/a = 400) the overlap between the 
fermions vanishes and the two peaks at E/t :^ 0.002 and 
E/t 0.007 merge into a single double-degenerate peak 
at energy E/t 0.001. The latter energy depends on the 
width of the nanowire^^. 



B. S-nanowire attached to a normal electrode 

In this Section we analyze the impact of a normal lead 
attached to the S-nanowire. This situation can arise for 
instance when the nanowire is only partially in proximity 
to a bulk superconductor so that part of the nanowire is 
in the normal state (as shown in Fig. 2). In order to 
get rid of finite-size effects, we consider a semi-infinite S- 
nanowire coupled to a normal lead. The case of a finite- 
length S-nanowire coupled to two normal leads at both 
ends does not yield additional significant information. 

In Fig. 5 we plot the LDOS, at different positions, 
as a function of energy for several values of the barrier 
strength 7. The two plots on the left [(a) and (c)] are for 
a non- trivial nanowire with /i/t = and V/t = 0.2, while 
the two plots on the right [(b) and (d)] refer to a topo- 
logically trivial nanowire with /i/t = and V/t = 0.6. 
Moreover, the top panels [(a) and (b)] refer to a position 
close to the interface in the S-nanowire, while the bottom 
panels [(c) and (d)] to a position close to the interface in 
the normal lead. When 7 is small the LDOS in the S- 
nanowire presents a finite gap Eg which is just a fraction 
of the superconducting pairing A owing to the presence 
of the Zeeman field. Namely, E^ 0.03 t = 0.3 A for 
the non-trivial nanowire and Eg 0.018 t = 0.18 A 
for the trivial case. Within the gap. Fig. 5(a) shows a 
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sharp Majorana peak at zero energy which broadens as 
7 increases^^'^^, and eventuahy disappears when 7^1. 
Such a peak can be fitted with the following Lorentzian 
function 



{T/2f 



(12) 



where r/2 is the half width at half maximum and M is 
the height at zero energy. The result of the fit is reported 
in Fig. 6: r/2 and M are plotted as functions of 7. Re- 
markably, r/2 depends only very weakly on the position 
in the S-nanowire where the LDOS is calculated. 

Interestingly, the Majorana peak is present also in the 
LDOS of the normal lead [Fig. 5(c)] as long as 7 is not 
exactly zero: the peak is still clearly distinguishable up 
to 7 :^ 0.6. Moreover, singularities develop at energies 
corresponding to ±£^g as 7 tends to 1. 

In the trivial phase, the LDOS in the S-nanowire at 
a position close to the interface [Fig. 5(b)] presents a 
single pair of peaks at 7 = (as compared to Fig. 4), 
since, being semi- infinite, the ends of the S-nanowire are 
sufficiently far away to be decoupled. Such peaks quickly 
broaden as 7 increases and eventually merge into a single 
peak at 7 0.3 (a further increase of 7 leads to the 
disappearing of the peak). As a result, as long as the 
coupling between the S-nanowire and the normal lead 
is not too strong, non-trivial and trivial phases can be 
distinguished from a measurement of the LDOS in the S- 
nanowire given a sufficiently large energy resolution (in 
the present case higher than 1% of the pairing A). On 
the contrary, only a very weak double-peak structure is 
visible in the LDOS in the normal lead at a position close 
to the interface and for small values of 7 [see Figs. 5(d)]. 

We mention that important clues on the topological 
phase of a S-nanowire coupled to a normal electrode can 
be also retrieved using a different transport setup with 
respect to the one depicted in Fig. 2. Indeed, instead of 
considering the current from the STM tip to the sample, 
it would in principle be possible to study directly trans- 
port through the NS junction present in the nanowire. 
Even though a theoretical analysis of this configuration 
is beyond the scope of the present work, in Appendix B 
we report on the low-bias conductance of the NS junc- 
tion. We find that, in the tunneling limit, the approx- 
imate quantization of the conductance can be exploited 
to identify the topological phase of the nanowire^^'^^'^^. 
We also remark that, even for transparent barriers, clear 
signatures of the presence of Majorana fermions can be 
extracted from the transport properties of the NS junc- 
tion provided that a quantum point contact is present 
close to the interface^^. 

Let us now analyze how the LDOS changes when the 
S-nanowire is driven through a topological phase transi- 
tion^^. In Fig. 7(a) we plot the topological invariant Q 
[calculated using Eq. (7)] and the number of open chan- 
nels (A^oc) as functions of the chemical potential fi for 
a fixed Zeeman field V = 0.3 t. The S-nanowire goes 
through a transition, from the non-trivial (Q = —1) to 
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FIG. 6: (Color online) Fitted values [through Eq. (12)] of the 
half width at half maximum r/2 (blue circles) and height at 
zero energy J\f (red triangles) of the peak in the LDOS at 
a position just inside the superconducting nanowire in the 
non- trivial phase. Solid lines are just guides for the eye. 



the trivial (Q = +1) phase, at /i 0.026 t. Interest- 
ingly, the number of channels increases from 1 to 2 at a 
much smaller value of the chemical potential fi 0.01 t 
such that the non-trivial phase persists even in the pres- 
ence of 2 open channels. This observation is in apparent 
contradiction with the intuitive picture adopted in the 
literature to explain the phase diagram of superconduct- 
ing nanowires, i.e. that a pair of Majorana fermions at 
the ends of the wire is associated with each open channel 
and that pairs of Majorana fermions on the same end can 
couple and form complex (Dirac) fermions. Accordingly, 
there should be a single isolated Majorana fermion at 
each end of the nanowire whenever the number of open 
channels is odd. On the other hand, this intuitive picture 
should be treated with care, simply because one concept 
(presence of Majorana fermions) is related to a supercon- 
ducting wire while the other (number of open channels) 
to a normal one. Indeed, already in Ref. 17 it was no- 
ticed that the system can be in the topologically trivial 
phase even when A/qc is odd. Here, we are observing the 
complementary situation in which the topological phase 
persists when A^oc is even. We believe that the underly- 
ing explanation is the same for both cases: the failure of 
the intuitive picture reported above to explain the whole 
phase diagram. As a consequence, we remark that the 
topological invariant is not necessarily in a one-to-one 
correspondence with the parity of the number of open 
channels. In Fig. 7(b) the LDOS as a function of energy 
is shown at a position inside the S-nanowire {x = 50 a, 
measured from the interface, and y = W/2 = 5 a) when 7 
is close to zero. Different curves correspond to different 
values of the chemical potential, with a vertical offset 
proportional to /i. The Majorana peak splits into two 
Dirac-fermion peaks at finite energy as the chemical po- 
tential moves through the phase transition. Besides, we 
also observe that the effective gap Eg initially decreases 
for increasing /i, then vanishes at the phase transition. 
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FIG. 7: (Color online) (a) Topological invariant Q (red solid 
line) and number of open transverse channels A^oc (blue dash- 
dotted line) as functions of the chemical potential /i for a sys- 
tem with V/t — 0.3. A phase transition occurs at /j/t ^ 0.026. 
(b) LDOS at a position inside a semi-infinite superconduct- 
ing nanowire for different values of the chemical potential fi 
keeping V/t = 0.3 fixed. Results have been offset vertically 
for clarity, with an increasing value of the chemical potential 
from bottom (/i/t = 0) to top (/i/t = 0.05). 



and thereafter increases again. Similar results (concern- 
ing the differential tunneling conductance at one end of 
a nanowire) have been recently presented in Ref. 18. 



C. SNS structure 

Let us now consider two semi-infinite S-nanowires con- 
nected through a normal nanowire (N-nanowire). For 
simplicity we assume transparent barriers at the inter- 
faces and we set 



A{i) = A{x/a) 







X < -L/2 
\x\ < L/2 
X > L/2 



(13) 



i.e. we allow for a finite phase difference Acp = cp^ — cpi, 
between the superconducting paring in the right and left 
S-nanowires, L being the length of the N-nanowire. In- 
dependently of the topological phase of the system, An- 
dreev bound states (ABSs) appear as sharp peaks in the 
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FIG. 8: LDOS as a function of energy at a given position 
close to the interface inside the left superconducting nanowire. 
Sharp peaks corresponding to Andreev bound states appear 
both within the effective gap Eg (~ 0.007 t) and inside the 
continuum. These results refer to the following set of param- 
eters: Acp = 2.0, fi/t = 0, V/t = 0.8, and L/a = 0. 



LDOS (see Fig. 8) and the ABS spectrum can thus be 
reconstructed from the energies at which such peaks oc- 
cur. As we shall discuss at length below, what differs be- 
tween topologically trivial and non-trivial phases is the 
parity of the number of zero-energy crossings in the ABS 
spectrum, which is related to the presence or absence of 
such a crossing at A(p = tt (protected by fermion par- 
ity). This result is in agreement with similar calculations 
performed for a strictly one-dimensional system^^'^^, for 
a two-dimensional SNS junction^^, and for a quantum 
spin-Hall insulator sandwiched between superconducting 
leads^'^. 

We start by considering the short-junction limit, L <C ^ 
(where oc ta/ Eg is the effective superconducting coher- 
ence length), in which we have a single ABS for each 
open channel. L = is a particular case in which the N- 
nanowire is absent and we have a step-like jump from 
to (fB. in the phase of a S-nanowire. In Fig. 9 we show, for 
L = 0, the ABS energies as a function of Acp for a topo- 
logically non-trivial [panel (a) for ja/t = and V/t = 0.2] 
and trivial [panel (b) for fi/t = 0.15 and V/t = 0.1] S- 
nanowires. In both cases we have a single ABS (with 
its opposite energy counterpart). The main difference 
between the two cases is that the number of zero-energy 
crossings is odd in the non-trivial phase while it is even in 
the trivial case, in agreement with Ref. 13. Furthermore, 
even though the overall periodicity of the ABS spectrum 
is 27r in both cases, the periodicity of a single branch of 
the ABS spectrum is Air only for the non-trivial situa- 
tion. In particular, the change of Acp by 27r at a fixed 
energy leads to the swapping of an ABS with its charge- 
conjugate state which has opposite fermion parity. In 
terms of Josephson current this leads to the fractional 
Josephson effect^ and can be interpreted either as a 47r- 
periodicity or as a two-valuedness of the Josephson cur- 
rent. For simplicity in the following we will address the 
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topologically non-trivial ABS spectrum as 47r-periodic. 

It is now interesting to compare the plot in Fig. 9(a) 
with the ABS spectrum relative to a short, one- 
dimensional, SNS Josephson junction with + Wy su- 
perconducting order parameter. An expression for the 
latter has been obtained, under the Andreev approxima- 
tion (which assumes an order parameter much smaller 
than the Fermi energy), in Ref. 38: 



(a) 



£;abs = ±Aabs VT cos (A(/?/2) 



(14) 



where T is the transmission probability of the N region 
and Aabs is an effective order parameter for the ABS. 
The solid line in Fig. 9(a) shows the result of a best fit, 
with respect to Aabs, once T is set to 1. The agreement 
between the numerical results (red empty circles) and 
the fit is quite good, despite the fact that the S-nanowire 
behaves effectively as a {px + ipy)-wdiYe superconductor 
only in the limit of large Zeeman fields, with respect to 
A (in the present case V = 0.2t and A = 0.1 t)^^. In the 
inset we show that, contrary to a non-topological SNS 
junction^^, the zero-energy crossing survives even when 
the barrier between the superconductors is not transpar- 
ent (7 = 0.1). In accordance with Eq. (14), the only 
effect of a transparency T < 1 is to reduce the band- 
width of the ABS spectrum preserving the crossing at 
Acp = TT. Indeed, with respect to the non-topological case 
discussed for instance in Ref. 40, here the ABSs have dif- 
ferent fermion parity, i. e. they correspond to states with 
an even or odd number of fermions, respectively. Since 
fermion parity is a conserved quantity in this case, the 
two ABSs can not couple (giving rise to the opening of 
a gap) even when the transparency of the system is low- 
ered. 

Let us now address the case of larger Zeeman fields, 
where, however, the situation is complicated by the fact 
that by increasing V also the number of open channels 
increases. In Fig. 10(a) the ABS spectrum, as a function 
of A(/?, is shown for a topologically non-trivial S-nanowire 
with /i = and V = 0.8 t, where three transverse open 
channels are allowed in the absence of superconducting 
pairing. The spectrum presents three ABS branches and 
an odd (three) number of crossings at zero energy, as 
expected for the non-trivial case. We fit separately the 
corresponding three ABS branches with Eq. (14) allowing 
three different values of Aabs, to be interpreted as effec- 
tive gaps for each open channel, one for each branch. As 
shown in the plot, the fits (solid lines) apparently approx- 
imate quite well the numerical data (red empty circles). 
The agreement, however, is completely lost for energies 
close to zero [see the inset of Fig. 10(a)]: i) only one 
ABS branch crosses zero energy at A(p = tt, ii) the spec- 
trum shows complicated avoided crossings which shift the 
zero-energy crossing of the two remaining ABS branches 
to A(f TT. As already discussed in Sec. IIIB, the Ma- 
jorana modes, which would localize at a given end if the 
three open channels were not coupled, strongly hybridize 
at Acp = TT giving rise to a true Majorana mode plus a 
pair of Dirac excitations at finite energy. Moreover, we 




7r/2 TT 37r/2 27r 
Acp 



(b) 




37r/2 27r 



FIG. 9: (a) Andreev bound state energies (red empty cir- 
cles) as a function of the phase difference Acp for topologically 
non-trivial superconducting nanowires (/i/t = 0, V/t = 0.2, 
and 7=1). The solid line is a best-fit using the expression 
in Eq. (14). Inset: same as in the main panel but for a non 
transparent barrier between the superconducting nanowires 
(7 = 0.1). (b) same as in the top panel but for topo- 
logically trivial superconducting nanowires {/n/t = 0.15 and 
V/t = 0.1). Both panels refer to a SNS junction with L/a = 0. 



notice that the periodicity of all three ABS branches is Att 
and the peaks in the LDOS connected to the two ABSs 
with larger energy appear (see also Fig. 8) even above 
the bulk gap £^abs (— 0.007 t in this case). 

The situation changes completely in the case of a 
smaller Zeeman field {V/t = 0.2), still in the topologi- 
cally non-trivial phase, while keeping the number of open 
channels equal to three {i.e. by setting ja/t = 1.7). The 
corresponding spectrum, reported in Fig. 10(b), presents 
again an odd (three) number of zero-energy crossings, 
but shows gap-like features: no fitting with Eq. (14) is 
now possible. This is due to the fact that now the ef- 
fective Px + ipy description no longer holds since ABSs 
with both chiralities are present^^, so that an inter-band 
s-wave pairing amplitude can mix them. In this case, one 
of the ABS branches has periodicity Air while the remain- 
ing two have periodicity 27r. We still expect a fractional 
Josephson effect. 

If we increase the length of the N-nanowire toward the 
long-junction limit {L ^ ^), the number of ABSs in- 
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7r/2 TT 37r/2 27r 



(b) 




0.03 



0.02 



0.01 



0.00 



-0.01 - 
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FIG. 10: (a) Andreev bound state energies (red empty circles) 
as a function of the phase difference /\Lp for topologically non- 
trivial superconducting nanowires when /j^/t = and V/t = 
0.8. The solid lines are fits to the function in Eq. (14) with 
T = 1 for three different values of Aabs for the three Andreev 
bound states. The inset shows a magnification of the region 
close to Acp = TT at small energies, (b) same as in the panel a 
but for /i/t = 1.7 and V/t = 0.2. Both panels refer to a SNS 
junction with L/a = 0. 



creases but the topological phase of the S-nanowire can 
still be detected from the number of zero-energy cross- 
ings in the interval Aip G [0, 27r]^^'^^. This is shown in 
Fig. 11 where the ABS spectrum is plotted for L/a = 50 
(long-junction limit) when the S-nanowires are in the 
non-trivial phase [panel (a) for V/t = 0.2 and /i/t = 0] 
and in the trivial phase [panel (b) for V/t = 0.05 and 
jii/t = 0]. As in the short-junction limit, the periodicity 
of at least one branch of the ABS spectrum is 47r only in 
the non-trivial case. 



IV. CONCLUSIONS 

In this paper we have considered a multi-band semicon- 
ducting nanowire subjected to spin-orbit coupling, super- 
conducting pairing and a longitudinal Zeeman field. De- 
pending on the values of such parameters, the nanowire 
presents a non-trivial topological phase in which a pair 
of Majorana modes, at an energy equal to the chemical 
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FIG. 11: (a) Andreev bound state energies (red empty circles) 
as a function of the phase difference Acp for topologically non- 
trivial superconducting nanowires (V/t = 0.2 and /n/t = 0). 
(b) same as in the top panel but for topologically trivial su- 
perconducting nanowires {V/t = 0.05 and /n/t = 0). Both 
panels refer to the long-junction limit {L/a = 50). 



potential, are localized at its ends. We have first derived 
an analytic expression for the phase boundaries of an in- 
finitely long multi-band nanowire. We have then numer- 
ically calculated and analyzed the local density-of-states 
of such nanowires in the case when they are coupled to 
normal regions (such as electrodes or links) and we have 
compared the topologically non-trivial and trivial phases 
in different situations. When the nanowire is coupled 
to a normal electrode we have found that the peak in 
the LDOS at zero energy (with respect to the chemical 
potential), corresponding to the Majorana mode, broad- 
ens with increasing coupling strength to the electrode, 
eventually disappearing for a transparent interface. In- 
terestingly, for finite coupling the peak is also present in 
the normal electrode, though being of smaller amplitude 
and broadening more rapidly with the strength of the 
coupling. In the trivial phase, and when the nanowire 
possesses two open channels in the absence of supercon- 
ducting pairing, a pair of peaks at finite energies appears 
as due to the hybridization of the two Majorana modes 
that would exist if the two channels of the nanowire were 
not coupled. Such peaks broaden with increasing cou- 
pling strength to the normal electrode, eventually merg- 



9 



ing for sufficiently large coupling. In the normal elec- 
trode only weak features survive. From the analysis of 
the topological phase transition, driven by varying the 
chemical potential at fixed Zeeman field, we have found 
that the nanowire remains in the topologically non-trivial 
phase even after the number of open channels goes from 
one to two. This suggests that, contrary to the intu- 
itive picture often referred to in the literature, the one- 
to-one correspondence between the topological invariant 
and the parity of the number of open channels is only 
approximate and should be treated with care. We have 
then considered the situation in which two semi-infinite 
nanowires (kept at different superconducting phases) are 
connected through a normal link of length L. Inde- 
pendently of the topological phase the density-of-states 
presents peaks due to Andreev bound states whose po- 
sition in energy depends on the superconducting phase 
difference /S.(p. While in the trivial phase the number 
of zero-energy crossings is even, in the topologically non- 
trivial phase this number is odd owing to the presence of a 
fermion-parity-protected crossing at /S.(p = tt. This differ- 
ence in the parity of the number of zero-energy crossings 
refiects the presence of at least one branch of Andreev- 
bound-state energy which is 47r-periodic (instead of the 
usual 27r-periodicity), leading to the so-called fractional 
Josephson effect. This anomalous 47r-periodicity of the 
Josephson current has been usually introduced in strictly 
one-dimensional systems while we have checked that it 
survives also in a multi-band nanowire, in agreement with 
Ref. 41. 
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Appendix A: Derivation of the phase diagram 

We first need to determine the eigenvalues and eigen- 
vectors of HBdcikx) for kx = and 7r/a, when V = A = 
0. The difference between kx = 0,7r/a is just a shift ±2t 
in the chemical potential, i.e. 

n^dGikx = 0, TT/a) = Ho + [so - (/i ± 2t)] r, . (Al) 

Moreover, we are assuming A = 0, so that 

^n^ 



Ho 



-(jynicry 



(A2) 



Thus, it suffices to consider the particle Hamiltonian Hp. 
For a wire of width W/a = n the characteristic polyno- 
mial Ln{£] Oi) of Hp can be defined recursively as 

Ln{e]a) = e Ln-i{e;a) - (t^ + Ln-2{£;a) , (A3) 



with Lo{s; a) = 1 and Li{e] a) = £. A formal solution to 
the recursive relation (A3) is given by 



with 



(A4) 



(A5) 



One can easily solve Ln{s] a) = and find the following 
expression for the energies ei 



cos 



Att 
n + 1 



A = l,...,n. (A6) 



For completeness, we mention that each eigenvalue is 
doubly degenerate and the corresponding eigenstates are 
characterized by an amplitude i/j^^ on site m and spin a 
given by 



^ e^^-^ sin (m ^ 
n + 1 



^^tl (A7) 



where the upper and lower signs refer to the two eigen- 
states and tan^ = a/t. Moreover, we observe that even 
allowing for a finite Zeeman field V 7^ 0, the eigenstates 
do not change since we are assuming the Zeeman field 
and the spin-orbit coupling in the transverse direction to 
be both proportional to . On the other hand, ifVy^O 
the two eigenstates in Eq. (A7) are no longer degenerate, 
but split to energies Si ± V. 

We now turn to the evaluation of the topological invari- 
ant Q in Eq. (7) when A 7^ 0. We need to compute the 
Pfaffians: 



Pf[HBdG(fex)^V^] with kx = 0,7r/a. 



(A8) 



We first introduce the 2n x 2n unitary matrix U of the 
eigenvectors (A7) of Hp and the associated 4n x 4n uni- 
tary matrix UbcIGi 



UBdG 



U 

iay W 



(A9) 



which diagonalizes HBdci^^^/d) when A = 0. Even if 
when A the following matrix 

^{kx) = Ul^^nBdG{kx)UBdG (AlO) 

is not diagonal, it is still useful to introduce it. Indeed, 
we can then write 

HBdG(fcx)<TV« = C/BdGl?(fcx)K^^C/BdG)^ (All) 

and, since by particle-hole symmetry we have that 

<TV^C/BdG = C/BdG^" , (A12) 
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the Pfaffian can be written as 



Pf[f/BdGl'(A:.)T"C/JdG] 

det(C/BdG)Pf[I>(fc.)r^ 
(-1)" Pf[P(fc,)T-] . (A13) 



The product inside the Pfaffian reads 



-iA U^aVU -E{k:,) 

lA WayU* E{k^) 
-E{k:c) -iA U^ayU 



1 

1 



(A14) 



where E(kx) is a diagonal matrix with entries given by 
{e\ + AV + £0 — M T A = ±1. This matrix is indeed 
antisymmetric. One can easily see that, upon a reorder- 
ing of the rows and columns described by a real unitary 
matrix V, the matrix V{kx)T^ can be put into a block 
diagonal form, where each block is a 4 x 4 matrix involv- 
ing the particle states with eigenvalues differing only by 
the sign of V and their hole counterparts. Namely, each 
block has the form 








-A 


A 







-s-V 





v 





-£ + 



■V 



V 



A 

/ 



(A15) 



where e = e\ + eq — ix^2t. Thus, this means that 

Pf[^BdG(fc.)<T^T''] = (-1)" Vi[V{k^)T-] (A16) 

= (-1)" det(V) X{[V^ -A^-{^i-ex-eo± 2tf] , 

A 

where the upper and lower signs refer to = and 
kx = 7r/a, respectively and the matrix V is the same for 
both kx = 0,7r/a. Thus, we finally have 



Q = sign{Pf[HBdG(0)aV^] Pf[H 

X,r]=±l 



BdG(^/a)^V^]} 
V 2tf - V^] 



(A17) 



and consequently the phase boundaries are given by 
Eq. (8). We finally mention that analogous results have 
been found in Ref. 42 in the case of spinless fermions in 
a p-wave superconducting nanowire. 



Appendix B: Transport across a NS junction 

In this Appendix we investigate another important tool 
to assess the topological phase of an S-nanowire coupled 
to a normal electrode: the low-bias conductance across 
a NS junction in the tunneling limit. In the limit of low 
bias, transmission through the superconductor is com- 
pletely suppressed and the conductance of the NS junc- 
tion can be expressed in terms of the Andreev reflection 
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FIG. 12: (Color online) Conductance (in units of 2e^ /h) of 
the NS junction as a function of the parameter 7 controlling 
the transparency of the barrier. The (blue) solid line refers to 
a S-nanowire in the topologically non-trivial phase {V/t — 0.2 
and 11 /t = 1.1), while the (red) dashed line to a S-nanowire 
in the topologically trivial phase {V/t = 0.2 and /x/t = 1.5). 



matrix (at the chemical potential) 



G 



2e^ 
h 



Tr 



2e^ 
h 



E 

m=l 



(Bl) 



Here Rm are the eigenvalues of the Hermitian matrix 
^Ite^he and A/qc is the number of open channels in the 
normal electrode. Owing to particle-hole symmetry, the 
Rm^ are either two- fold degenerate or equal to or 1^3,43 
The presence of a fully Andreev-reflected mode (giving a 
quantized contribution to the conductance) is a signature 
of the existence of an uncoupled Majorana fermion at the 
Fermi energy^^'^^. As a consequence it is possible to write 
the conductance in the following form^^ 



G 




Q^^Y^Rn 



(B2) 



where the primed sum is restricted to the degenerate An- 
dreev reflection eigenvalues and Q is the topological in- 
variant in Eq. (7). In the limit of poorly transparent 
barriers (7 <C 1), we expect that almost all modes are 
fully reflected {Rm ~ 0, though never exactly zero^^'^^) 
and thus 



2p2 



(B3) 



As a consequence the low-bias conductance of the NS 
junction in the tunneling limit gives an important infor- 
mation on the topological phase of the S-nanowire^^'^^'^^. 
This result can be extended to almost transparent barri- 
ers if we include a ballistic quantum point contact close to 
the NS interface^^. In Fig. 12 we show the conductance 
G (in units of 2e^ /h) as a function the parameter 7 con- 
trolling the transparency of the barrier. In the tunneling 
limit (7 <C 1), we notice that the conductance approaches 
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for a topologically trivial S-nanowire (red dashed line, non-trivial S-nanowire (blue solid line, V/t = 0.2 and 
V/t = 0.2 and /i/t = 1.5) or 2e^ /h for a topologically /i/t = 1.1), in agreement with Eq. (B3). 
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